CERN-TH-2002-023 



Critical exponents from optimised renormalisation group flows 

Daniel F. Litim 

Theory Division, CERN, CH - 1211 Geneva 23. 
(Daniel. Litiin@cern.ch) 

Abstract 

Within the exact renormahsation group, the scahng solutions for O(A^) symmetric 
scalar field theories are studied to leading order in the derivative expansion. The 
Gaussian fixed point is examined for d > 2 dimensions and arbitrary infrared reg- 
ularisation. The Wilson-Fisher fixed point in d = 3 is studied using an optimised 
flow. We compute critical exponents and subleading corrections-to-scaling to high 
accuracy from the eigenvalues of the stability matrix at criticality for all N . We es- 
tablish that the optimisation is responsible for the rapid convergence of the flow and 
polynomial truncations thereof. The scheme dependence of the leading critical ex- 
ponent is analysed. For all > 0, it is found that the leading exponent is bounded. 
The upper boundary is achieved for a Callan-Symanzik flow and corresponds, for all 
N , to the large- limit. The lower boundary is achieved by the optimised flow and 
is closest to the physical value. We show the reliability of polynomial approxima- 
tions, even to low orders, if they are accompanied by an appropriate choice for the 
regulator. Possible applications to other theories are outlined. 



I. INTRODUCTION 



Renormalisation group techniques are important tools to describe how classical physics 
is modified by quantum fiuctuations. Integrating-out all quantum fiuctuations provides the 
link between the classical theory and the full quantum effective theory. Universality implies 
that the details of the underlying classical theory - other than the global symmetries, long- 
or short-range interactions, and the dimensionality - are irrelevant for the characteristics of 
the quantum effective theory. For this reason, universal properties of phase transitions in 
numerous physical systems (entangled polymers, liquid-vapour transition, superfiuid tran- 
sition in ^He, ferromagnetic transitions, QCD phase transition with two massless quark 
fiavours) can be addressed based on simple scalar field theories [||]. 

A useful method is given by the Exact Renormalisation Group (ERG) [0-0, which is 
based on the Wilsonian idea of integrating-out infinitesimal momentum shells. The corre- 
sponding fiow, which has a simple one-loop structure, is very fiexible concerning approxi- 
mations, and its domain of applicability is not tied to weak coupling. Recently, it has been 
shown that ERG fiows can be optimised, thereby providing improved results already to low 
orders within a given approximation ||?|-p!0|. In the present paper, we apply this idea to the 
universality class of 0{N) symmetric scalar theories in three dimensions and compute crit- 
ical exponents and subleading corrections to scaling. We expect that insights gained from 



this investigation will also prove useful for applications to gauge theories [TTI] or gravity 0] , 
which are more difficult to handle. 

Universal critical exponents have been computed previously using either polynomial 
truncations of exact renormalisation group fiows, or the derivative expansion to leading 



and subleading order ||10| , p!3| -p^ , and in p^ , |25[| based on the proper-time renormalisation 
group All results are affected by the underlying approximations which induce a spu- 
rious dependence on the regularisation pi]|jr7| , |2^P7| -pi| . This is somewhat similar to the 
scheme dependence within perturbative QCD, or within truncated solutions of Schwinger- 
Dyson equations. While this scheme dependence should vanish at sufficiently high order 
in the expansion, practical applications are always bound to a finite order, and hence to 
a non-vanishing scheme dependence. In some cases, it has even been observed that higher 



order results happen to be worse than lower order ones In consequence, one should 



gain some understanding of the spurious scheme dependence. Without this, it is difficult to 
decide which of the different scheme-dependent results within a fixed truncation could be 
considered as trustworthy. 

A partial understanding of the interplay of approximations and scheme dependence has 
been achieved previously. For scalar QED the scheme dependence in the region of 



first order phase transition has been studied in 123123]. For 3d scalar theories, the interplay 



between the smoothness of the regulator and the resulting critical exponents has been ad- 



dressed in using a minimum sensitivity condition. For Einstein quantum gravity, where 
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a new UV fixed point has been found recently to low orders in a polynomial truncation, the 
corresponding analysis has been given in |^|. The weak scheme dependence found in these 



cases suggests that higher order corrections remain small, thereby strengthening the results 
existing so far. The evidence created this way is partly circumstantial, because the range 
over which a physical observable varies with the scheme depends on the class of regulators. 

Here, we address the problem from a different perspective. The main ingredient in our 
analysis is the concept of optimisation HTHToj. For a given physical problem, it should be 



possible to identify specific regulators which lead to a better convergence behaviour of the 
flow. This strategy is based only on the ERG flow itself and has lead to a simple opti- 
misation criterion for flows 0. Optimised flows have a number of interesting properties 
PJ. They lead to a fast decoupling of heavy modes, they disentangle quantum and thermal 
fluctuations along the flow, and they lead to a smooth approach towards a convex effective 
potential for theories in a phase with spontaneous symmetry breaking. The optimisation is 
closely linked to a minimum sensitivity condition |^ in a sense which will be made transpar- 
ent below. Furthermore, optimised flows have been shown to improve the convergence of the 
derivative expansion |10|. Thus, optimised flows are promising candidates for high precision 



computations within this formalism. Here, we do so within a local potential approximation. 

The second new ingredient of our analysis consists in a study of the largest possible 
range of flows, and the corresponding critical exponents. We find that the range is larger 
than previously assumed. Furthermore, the results from optimised flows are located at the 
(lower) boundary and happen to be closest to the physical values. 

For the numerical analysis, and apart from the local potential approximation, we em- 
ploy a polynomial approximation. This additional approximation is reliable if it converges 
reasonably fast towards the full solution. However, it has been criticised previously in the 
literature. For a sharp cut-off, it has lead to spurious solutions , and its convergence was 
found to be poor |3^, which has lead to strong doubts concerning its reliability (see also 
|33|). In contrast to these findings, we show that the poor convergence is an artifact of the 



sharp cut-off regulator, rather than an artifact of the polynomial approximation. Using an 
optimised flow, we find that the polynomial approximation is stable and that it converges 
rapidly for all technical purposes. 

The format of the paper is as follows. We review the basic ingredients of the formalism 
and introduce the optimisation ideas (Sect. |ID. Then, we introduce our numerical method 
and study the non-trivial scaling solution in 3d (Sect. We compute the eigenvalues at 
criticality to high accuracy from an optimised flow. The convergence and stability of the 
flow, and of the polynomial truncation, are established (Sect. [TVI ). We study the scheme 
dependence of the critical index u (Sect. |V|). Finally, we discuss the main results of the 
paper with particular emphasis on the predictive power, on the convergence properties of 
flows, and on implications for other theories (Sect. [V^) . Three appendices contain the study 
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of the Gaussian fixed point for arbitrary regularisation and d > 2 (App. and technical 
details for specific classes of regulators, including a computation of the corresponding flows 
and critical exponents (Apps. ^ and y). 

II. RG FLOW FOR 0{N) SYMMETRIC SCALAR THEORIES 

Here, we briefly review some basic ingredients of the ERG formalism, and its approxi- 
mation to leading order in the derivative expansion. We also discuss important aspects of 
the regularisation, and its optimisation, which is employed in the following sections. 



A. Renormalisation group flows 

Exact renormalisation group equations are based on the Wilsonian idea of integrating 
out momentum modes within a path integral representation of quantum field theory [Q. In 
its modern form, the ERG flow for an effective action for bosonic fields is given by the 
simple one-loop expression 




dtTk[(i)] = -Tr — + Rk dtRk (2.1) 



Here, t = In is the logarithmic scale parameter, and Rk{q^) is an infrared (IR) regulator at 
the momentum scale k. From now on, we suppress the index k on R. The flow trajectory of 
Eq. (|2.1| ) in the space of action functional depends on the IR regulator function R. R obeys 
a few restrictions, which ensure that the flow equation is well-defined, thereby interpolating 
between an initial action in the UV and the full quantum effective action in the IR. We 
require that 

hm i?(g2)>0, (2.2) 
hm i?(g2)^0, (2.3) 
lim R{q^) oo . (2.4) 

fc— >A 

Equation ( |2.2| ) ensures that the effective propagator at vanishing field remains finite in the 
infrared limit — 0, and no infrared divergences are encountered in the presence of massless 



modes. Equation (|2.3| ) guarantees that the regulator function is removed in the physical 
limit, where F = limfc^oTfe- Equation (|2.4|) ensures that F^ approaches the microscopic 
action S = limk^^Tk in the UV limit k —>■ A. We put A = oo in the sequel. For later use, 
we introduce a dimensionless regulator function r{y) as 

R{q^)=q^r{q^/k^). (2.5) 
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Now we turn to the flow equation for an 0{N) symmetric scalar field theory in d dimensions 
to leading order in the derivative expansion, the so-called local potential approximation [P^] . 
This approximation amounts to the Ansatz 

r, = J d'^x (UuiP) + \d,rd^,(i>^ . (2.6) 

for the functional F^. The Ansatz neglects higher order corrections proportional to the 
anomalous dimension t] of the fields. The latter are of the order of a few percent for the 
physically interesting universality classes = 0, ■ ■ ■ , 4. Hence, we expect that a derivative 
expansion is sensible, and that the result of a leading order computation is correct up to 



corrections of the order of rj. Inserting this Ansatz into (p.l|) and evaluating it for constant 
fields leads to the flow for Uk- We rewrite this flow equation in dimensionless variables 
u{p) = Uk/k'^ and p = In addition, the angular integration of the momentum 



trace is performed to give 



dtu + du-{d- 2)pu = 2vd{N - l)i{u') + 2vdi{2pu") (2.7) 
with v^^ = 2'^+-'^7r'^/^r(|). The function i{uj) are given by 

^ ^ ^Jo y{l + r)+uj ^ ^ 



with y = q^/k"^ and dtr[y) = —2yr'[y). The flow ( p.7|) is a second order non-linear partial 
differential equation. All non-trivial information regarding the renormalisation flow and the 
regular isat ion scheme (RS) are encoded in the function ( p.8|) . The momentum integration is 
peaked and regularised: for large momenta due to the regulator term dtr{y), and for small 
momenta due to r(?/) in the numerator. All terms on the left-hand side of Eq. ( |2.7] ) do not 
depend explicitly on the RS. They simply display the intrinsic scaling of the variables which 
we have chosen for our parametrisation of the flow. 



B. Optimisation 



A good choice for the regulator is most important for a rapid convergence and the stabil- 
ity of an approximated flow towards the physical theory. Recently, is has been argued that 
such (optimised) choices of the IR regularisation are indeed available J^^- The main ob- 
servation is that the flow trajectory of ( p.l| ) depends on the regularisation. This observation 
is most important for approximated flows: typically, their endpoint also depends spuriously 
on the regularisation. The dependence is absent for the full integrated flow. Hence, in order 
to provide reliable physical predictions, it is important to seek for regularisations for which 
the main physical informations are already contained within a few leading order terms of 
an approximation. This issue is intimately linked to the stability of flows. 
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The main ingredient in is the full inverse propagator. Due to the IR regularisation, 
the full inverse regularised propagator displays a gap as a function of momenta MM , 



mm 

g2>0 



[64>{q)64>{-q) 



+ Rk{q') 1 =CF>0. (2.9) 



The functional derivative is evaluated at a properly chosen expansion point 0o- The ex- 
istence of the gap C > implies an IR regularisation, and is a prerequisite for the ERG 
formalism. Elsewise, ( p^.l] ) becomes singular at points where the full inverse effective prop- 
agator develops zero modes. It is expected that an approximated ERG flow in the space of 
all action functionals is most stable if the regularised full propagator is most regular along 
the flow. This reasoning corresponds to maximising the gap. 

To leading order in the derivative expansion, and dropping irrelevant momentum- 
independent terms, the gap is given by 

C = miny(l + r(2/)). (2.10) 

Within the local potential approximation, the gap C in ( p.lOl ) depends on the regularisation, 
but not on the specific theory considered. A comparison of the gap of different regulators 
requires an appropriate normalisation of r. In order to make the flow (|2.1|) more stable, 
we require that the g 5 clS Cl function of the RS, is maximal in the momentum regime 
where the flow receives its main contributions. This is the optimisation criterion of 
It corresponds to an optimisation of the radius of convergence of amplitude expansions and 
the derivative expansion. It can also be shown that the optimisation is closely linked to a 
minimum sensitivity condition. Optimised regulators are those for which the maximum in 
( [^.1(J| ) is attained. In Ref. these considerations have lead to a very specific solution to 
the optimisation condition, given by 

i?opt(g') = (A:'-g')^(fc'-g'). (2.11) 

It has a number of interesting properties. For large momenta > k'^, the propagation of 
modes is fully suppressed since R = 0. For small momenta q^ < k"^, all modes propagate 
with an effective mass term given by the IR scale, + R{q^) = k"^ ■ Based on the discussion 
in Refs. we expect that Eq. (|2.11|) leads to an improved convergence and hence better 



physical predictions already to leading order in the derivative expansion. Rewriting ( p.ll| ) 
in the form of ( |2.5| ) leads to 

ropt(y) = (i-l)^(l-y) (2.12) 

Below, we mainly employ this regulator, and variants thereof (cf. Appendices and 0). 
When expressed in terms of the optimised regulator (|2.11|) , the fiow equation (|2.7| ) becomes 
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4 - 1 4 1 
dtu = -du + {d- 2)pu' + -;Vd-—^ + -:v. ^ „ . (2.13) 

Notice that the numerical factors (4t>rf) / d can be absorbed into the potential and the fields 
by an appropriate rescaling. 



III. FIXED POINTS 



The fiow equation (|2.13|) is known to exhibit two scaling solutions in 2 < c/ < 4, which 
correspond to the Gaussian and the Wilson-Fisher fixed point, respectively. The Gaussian 
fixed point is given by the trivial solution u^, = const, and is discussed in the Appendix ^ 
for arbitrary regularisation and 2 < d < 4. In this section, we study the non-trivial fixed 
point 7^ const, in Zd based on an optimised fiow. We introduce the numerical method, 
and discuss the scaling form of the fixed point solution. Universal critical exponents are 
computed in the next section. 



A. Numerics 

Numerical methods for solving partial differential equations are well-known. Here, we 
employ a polynomial truncation of the scaling potential, retaining vertex functions 0^" up 
to a maximum number ntrunc- Polynomial approximations are reliable if they depend only 
very weakly on higher order operators beyond some finite order of the truncation. (In 
the following section, we explicitly confirm that this is indeed the case.) Two different 
polynomial expansions of the potential are used: expansion I corresponds to 

"^trunc 1 

«(P)= E -^nP"-. (3.1) 

In Eq. ( p.l|) , Ai denotes the (dimensionless) mass term at the origin. We have normalised the 
potential as m(p = 0) = 0. All higher order coefficients A„ denote the n-th order coupling at 
vanishing field. Expansion II, alternatively, approximates the potential about the potential 
minimum p = po as 

ntrunc 1 

«(P)= E ^A„(p-Ai)". (3.2) 

In Eq. (|3.2| ), Ai denotes the location of the potential minimum defined by u'{p = Ai) = 0. 
We have normalised the potential as u{Xi) = 0. All higher order coefficients A„ denote 
the n-th order coupling at the potential minimum. Both expansions ( p.l] ) and ( p.2| ) are 
symmetric under —>■ —0, and approximate the potential as an even polynomial in 0. 
The number of independent operators contained in Eq. p.l| ) or Eq. ( p.2| ) is ntrunc- The 
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expansions transform the partial differential equation ( |2.7| ) into ritrunc coupled ordinary 
differential equations dtXi = j3i{{Xn}) for the set of couplings 

{A„, 1 <n< ritrunc} ■ (3.3) 

For the numerical study, we use mostly expansion II, which is known to have superior 
convergence properties in comparison to expansion I |]33|f?[|. 



B. Wilson-Fisher fixed point 



The Wilson-Fisher fixed point corresponds to the non-trivial scaling solution of (|2.13| ). 



Here, we restrict ourselves to the case d = 3, and to the optimised regulator as introduced 
above. The scaling potential obeys the differential equation 

1 iV- 1 1 1 , , 

= — 3m* + , H 1 h const. (3.4) 

^ * 67r2 1 + < 67r2 1 + < + 2pu'l ^ ' 



and u*(p) 7^ const. The constant in Eq. (|3.4|) can be fixed freely, and we chose it such 
that the scaling potential vanishes at its minimum, -u* = at the point po with u'^ = 0. 
An analytical solution of Eq. (|3.4|) has been given in Ref. [T^ for the limit N = oo. For 



7^ oo, and in the vicinity of p = 0, the scaling solution can be obtained analytically as a 
Taylor expansion in the field. For d = 3 and = 1, and inserting the expansion ( p.l|) into 
]4D with const. = 0, we find 

Ao= j^7r-2(l + Ai)-i (3.5a) 

A2 = -47r2Ai(l + Ai)^ (3.5b) 

A3= fivr^Ai(l + Ai)' (l + 13Ai) (3.5c) 

A4 = -IfSyr^ A? (1 + Ai)' (1 + 7 Ai) (3.5d) 

A5= ^7r«A?(l + Ai)'(2 + 121Ai + 623A^) (3.5e) 



for small fields. Similar explicit solutions are found for 7^ 3 or A^ 7^ 1. Notice that all 
couplings are expressed as functions of the dimensionless mass term at vanishing field, Ai. 
The domain of validity of this expansion is restricted to < p < pc < 00. From the explicit 
solution for large fields, ( p.6|) below, it is clear that the polynomial approximation cannot 
be extended to arbitrary large fields. The value pc defines the radius of convergence for 
the polynomial approximation. It is linked to the gap parameter C introduced earlier [0]. 
Not all values for Ai lead to a scaling solution which remains finite and analytical for all 
p < Pc- It is at this point where the quantisation of Ai becomes manifest: only two values for 
Ai correspond to well-defined solutions of the fixed point equation, given by the Gaussian 
fixed point Ai = and the Wilson-Fisher fixed point Ai = Ai* < 0. The value for Ai* is 
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determined by fine tuning Ai such that the solution ( |3.5[ ) extends to p = pc- The result is 
Ai^ = —0.1860642 ■ ■ ■ for = 1. In the other limit, for large fields p ^ 1, we find 

u{p)=Ap^ + ^A~'p~' + ..., (3.6) 

(and similarly for 7^ 1), where the dots denote subleading terms in p and the constant 
A > has to be fixed appropriately, in a way similar to Ai in ( p.5|) . It would be interesting 
to study the analyticity properties of the fixed point solutions more deeply, and to contrast 
it with the analysis of |^ for the sharp cut-off flow. Despite the simple explicit form of the 
flow, and its explicit solution for large and small fields, it is more efficient to identify the 
scaling solution and the related critical exponents using numerical methods. 

C. Scaling potential 

We have solved the resulting coupled set of differential equations in c? = 3 dimensions, 
for = 0,1, ■■■,10 and for truncations up to ntrunc = 20. We have retained only the 
solution which has one unstable eigendirection, corresponding to the Wilson-Fisher fixed 
point. Fig. 1 shows the corresponding scaling potential Ui, in the vicinity of the potential 
minimum. In Fig. 1, the unusual normalisation of the potential has been chosen only for 
display purposes. The scaling potential has, for all A^, a minimum at non-vanishing field. 
For this reason, the expansion II has better convergence properties than expansion I. 




Figure 1: The scaling potential, with cat = Figure 2: The amplitude u'^(/9) of the scaling 
40 - 2N. From left to right: = 0, 1, • • • , 10. solution. Prom left to right: A^ = 0, 1, • • • , 10. 
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Fig. 2 shows the amphtude u' at the fixed point. It is a monotonic function of the field. 
It approaches its most negative values at vanishing field. The fiow equation would have 
a pole at points where 1 -\- u' vanishes. This is, however, never the case, because 1 -\- u' 
stays always positive. Fig. 3 shows the radial mode of the scaling potential u'^ + 2f>u'l in 
the vicinity of the potential minimum. It is a monotonic function of the fields. The dots 
in Fig. 3 indicate where the derivative u'^ changes sign. Again, the most negative value is 
attained at vanishing field and decreases with increasing iV, but it stays always above the 
pole of the fiow equation, l-\-u'^ + 2pM^' > 0. 

The scaling solution is non- universal. However, critical exponents or the eigenvalues of 
small perturbations about the scaling solutions are universal. For their determination, it is 
sufficient to study the fiow of small perturbations in the vicinity of the scaling potential, 
which is done next. 

IV. CRITICAL EXPONENTS 

In this section, we compute the critical exponents to leading order in the derivative ex- 
pansion. Numerical results are given up to six significant figures (Tab. 1). A higher precision 
can be achieved, but is not required at the present state. We find a rapid convergence of the 
polynomial approximation, for all observables considered. Our results are compared with 
those from other regulators. 
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A. Eigenvalues 



Given the Wilson-Fisher fixed point solution, we can seek for universal critical exponents. 
In the vicinity of the non-trivial fixed point, we have to solve the eigenvalue equation 



(4.1) 



in order to determine the various universal eigenvalues u. Using the fiow equation, setting 
d = 3 and choosing m = 1, and expanding the fiow to leading order about the Wilson-Fisher 
fixed point, we find 







+ 



LU + 2- 

1 



67r2 (1 + u'^y 
P 



6u' + 



37r2 



2N 
l + < 



l+<-p<-2pV" 
(l + < + 2p<)2 



6u" 



37r2 1 + < + 2fm'l 



6u ' 



(4.2) 



Instead of solving Eq. ( ^.2| ) directly for an eigenperturbation 6u' with eigenvalue u, we follow 
a slightly different line which is numerically less demanding. Based on the polynomial 
expansion used in the previous section, the fixed point potential is given by the set of 
couplings {A„^^}. At the fixed point, the fiow of the couplings dfXn = Pn vanishes, Pn{K,-k) = 
0. The eigenvalues u of the stability matrix at criticality 



Mi, 



J Ia=a* 



(4.3) 



correspond to the eigenvalues of Eq. (|4.1| ). Hence, the problem of finding the eigenvalues of 
Eq. ([4.1|) reduces to the problem of finding the eigenvalues of the stability matrix M. 



N 


V 




<^2 









0.592083 


0.65788 


3.308 


6.16 


9.2 


1 


0.649562 


0.655746 


3.180 


5.912 


8.80 


2 


0.708211 


0.671221 


3.0714 


5.679 


8.440 


3 


0.761123 


0.699837 


2.9914 


5.482 


8.125 


4 


0.804348 


0.733753 


2.9399 


5.330 


7.867 


5 


0.837741 


0.766735 


2.9108 


5.2195 


7.665 


6 


0.863076 


0.795815 


2.8967 


5.1409 


7.512 


7 


0.882389 


0.820316 


2.8916 


5.0863 


7.396 


8 


0.897338 


0.840612 


2.89163 


5.04848 


7.3086 


9 


0.909128 


0.857384 


2.89438 


5.02232 


7.2425 


10 


0.918605 


0.871311 


2.89846 


5.00420 


7.1921 


oo 


1 


1 


3 


5 


7 



Table 1: The first five eigenvalues < < wi < 0^2 < '^s < "^4, with v 
N and d = 3 dimensions. 



-1/iOo > for various 
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B. Results 



We have computed the eigenvalues of (|4.3| ) as functions of the order of the polynomial 
truncation. The numerical results for the first five eigenvalues for various N and d = 3 
dimensions are given in Tab. 1 and in Figs. 4-8. Results for other regulators are compared 
in |jlO[ (see also Sects. [IV D| and [VD , and results for the asymmetric corrections to scaling are 
given elsewhere. Exact results, independent on the regularisation, are known for N = oo. 
The exact eigenvalues are given by Un = 2n—l+0{l/N), and only the subleading corrections 
0{1/N) depend on the regulator. For N = —2, it is known that u = ^. Tab. 1 shows the 
first five eigenvalues, ordered as uq < < u < uj2 < uj^ < uj^. The critical exponent v is 
given by z/ = — I/cuq > 0. As a function of A^, the eigenvalue v interpolates monotonically 
between the exact values v = \ iox N = —2 and z/ = 1 for = oo. The eigenvalues and 

are non-monotonic functions of A^. For small A^, lo decreases until A^ 1, and increases 
towards the exact asymptotic value a; = 1 for N ^ oo. The eigenvalue 002 decreases until 
A^ ^ 7 — 8 before settling to the asymptotic value C1J2 = 3 at A^ = 00. The eigenvalues uj^ 
and LJ4 are monotonically decreasing functions of A^. 



0.9 



0.8 



0.7 



0.6 




0.5 

4 6 8 10 12 14 16 18 

^trunc 

Figure 4: The exponent ^{N) as a func- 
tion of N and of the order of the truncation. 
From top to bottom: N = 10, 4, 3, 2, 1, 0. 



0.7 



UJ{N) 



0.6 \ 

6 8 10 12 14 16 18 20 

^trunc 

Figure 5: The eigenvalue uj{N) func- 
tion of N and of the order of the truncation. 
From top to bottom: N = 10, 4, 3, 2, 0, 1. 



Fig. 4 shows the results for z/ as a function of ntrunc and A^. It is seen that the expan- 
sion is very stable. It convergences already within low order of the truncation towards the 
asymptotic value. Typically, ritrunc ~ 6 gives the correct result below the percent level. Fur- 
thermore, the convergence is better for larger values of A^. This behaviour is observed for all 
eigenvalues (e.g. Figs. 4-8). Fig. 5 shows the same behaviour for the smallest subleading 
eigenvalue u. For ritrunc ~ 8, it has settled below the percent level of the correct result. 
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3.6 



3.4 



3.2 



2.8 



2.6 




U2{N) 



8 



14 16 18 



10 12 
^trunc 

Figure 6: The critical index UJ2{N) as a func- 
tion of N and of the order of the truncation. 
Prom top to bottom: N = 0,1, 2, 3, 4, 10. 



6 8 10 12 14 16 18 20 
^trunc 

Figure 7: The critical index U!3{N) as a func- 
tion of N and of the order of the truncation. 
Prom top to bottom: N = 0,1, 2,3, A, 10. 




8 10 12 14 16 18 20 
^trunc 



Figure 8: The critical index UJ4{N) as a func- 
tion of N and of the order of the truncation. 
Prom top to bottom: iV = 0, 1, 2, 3, 4, 10. The 
isolated point at n = 7 corresponds to = 10, 
and the two at n = 9 to A?^ = 2 (upper) and 
N = 3 (lower). The intermediate points (at 
n = 8 and n = 10, resp.) are missing because 
they have a small imaginary part. 



For low orders of the truncation, the eigenvalues U2, UI3 and a'4 (Figs. 6, 7 and 8, respec- 
tively) depend more strongly on ritrunc, and do not even lead to real eigenvalues in some cases 
(c.f. Figs. 7 and 8). Again, the dependence is even stronger for smaller N. For sufficiently 
high order in the truncation, however, all eigenvalues are real, and the convergence towards 
the asymptotic value is fast. Typically, U2, co^ and 004 reach their asymptotic values below 
the percent level for ritrunc ~ 10, 12 and 14. 

For all eigenvalues, the basic picture is the same: for small ritruno the truncation tends 
to overshoot the asymptotic value, but with increasing ntrunc it relaxes towards it with a 
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remaining oscillation and decreasing amplitude. From a specific order onwards, the trun- 
cation sits - for all technical purposes - on top of the asymptotic value. For a fixed level 
of accuracy, a lower order in the truncation is required for the dominant observables like 
u or uj, while a higher order is required for the sub leading eigenvalues. Roughly speaking, 
to obtain the eigenvalue uJn,n = 0, ■ ■ ■ accurate below the percent level, a truncation with 
'^trunc ~ 2n + 6 independent couplings is required. 



C. Convergence and stability 



Next, we discuss the convergence and stability of the polynomial approximation for 
an optimised flow. From the results presented so far (cf. Figs. 4 - 8), we conclude that 
the optimised flow ( |2.13[ ) leads to a fast convergence of the polynomial approximation for 
the scaling potential. More importantly, we have seen that the inclusion of further vertex 
functions — increasing ritrunc — ^ ^trunc + 1 — does not alter the fixed point structure. Rather, 
it leads to a small modification of the actual fixed point solution and to minor corrections 
for the critical exponents. This implies that the flow is very stable, and that most of the 
physical information is already contained in a few leading order terms of the truncation. 
This result is by no means trivial. A counter example is furnished by the sharp cutoff (see 
also the following section), where the convergence of the polynomial approximation is poor. 
Here, the good convergence hinges on the use of an appropriately optimised regulator 




Figure 9: Ising universality class. Conver- 
gence of ftrunc (expansion II) towards fopt with 
increasing truncation. Points where ftrunc is 
larger (smaller) than z/opt are denoted by o (•). 
Roughly speaking, for 2 < ntrunc < 20, the ac- 
curacy of the critical exponents improves by 
one decimal point every An 2 — 2.5. 



10 15 

^trunc 



Let us have a closer look at the rate of convergence towards the asymptotic values of 
expansion II (cf. Fig. 9). We denote with Utmnc the approximate critical exponent which 
retains n = ntrunc independent parameters in the effective action. The semi-logarithmic 
plot in Fig. 9 then shows the rate with which successive approximations converge towards 
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the asymptotic value Uopt- The series z/trunc oscillates about the asymptotical values with 

a decreasing amplitude and, roughly, a four-fold periodicity in the pattern + H . The 

curve in Fig. 9 can be approximated by a straight line with a slope ~ —.4 to —.5. Hence, 
for every An ^ 2 — 2.5, the accuracy of i^trunc increases by one decimal place. Here, we have 
analysed the convergence for = 1. For larger A^, the convergence is typically faster than 
for N = 1, while for = 0, it is about the same. Hence, the present considerations are 
qualitatively the same for all A^. 
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Figure 10: The critical index v for the Ising universahty class as a function of the order 
of the truncation ntruno for expansion I around vanishing field and expansion II around 
the non-trivial minimum. Left panel: I and II converge towards the asymptotic value. 
Right panel: Magnification by a factor of 625. The variation of all data points is below 
10^'^. II converges faster than I. Expansion I (II) fluctuates about the asymptotic value 
with decreasing amplitude and a six-fold (four-fold) periodicity; see also Fig. 9. 



Until now, we have only employed expansion II, defined in (|3.2|) . It is worthwhile to 
employ as well expansion I, defined in ( |3.1| ). This has been done for the critical exponent u 
of the Ising universality class in Fig. 10. The left panel shows that both expansions lead to 
a fast convergence towards the asymptotic value. The right panel is a magnification by 625, 
showing that expansion II indeed converges faster, although for ritrunc = 10 their difference 
is already below 10"'^. Expansion I fluctuates about the asymptotic value with decreasing 
amplitude and a six-fold periodicity in the pattern -|- -f- H . 

D. Convergence and scheme dependence 

In this section, we discuss the convergence of the ERG flows and the polynomial expan- 
sion for various regulators. In Fig. 11, we have computed the critical exponent v [N = 1) 
for the sharp cutoff Vsharpiy) = 1/ ^(1 — u) — ^, the quartic regulator rquart(y) = y^^ and the 
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optimised regulator ( |2.12| ). The left (right) panel uses the expansion I (II). Both rquart and 
Topt are optimised regulators in the sense coined in section IIB . From Fig. 11, three results 
are noteworthy. First, for the expansion I, we confirm that the convergence is very poor for 
the sharp cutoff. For both the quartic and the optimised regulator we find a good conver- 
gence. Second, the convergence is additionally improved by switching to the expansion II. 
Third, the critical exponents obey z/sharp > t'quart > t'opt- Hence, the better the convergence 
and the stability of the flow, the smaller the resulting critical exponent u. This observation 
is also linked to the convergence of the derivative expansion [[10 . 




Figure 11: The critical index u for the Ising universality class. Results are given for 
Expansion I (left panel) and II (right panel), and for the sharp cutoff (upper curves) the 
quartic regulator (middle curves) and the optimised regulator (lower curves). 

In 0, it has been shown that the gap C (the radius of convergence for amplitude ex- 
pansions) is linked to the radius of convergence C for the expansions ( p.l| ) and (|3.2| ). An 
alternative way for identifying C consists in studying the complex structure of the scaling 
solution. For real p, the scaling solution is finite and real for all p < oo. In the complex 
plane, the scaling solution has (various) poles. Those closest to the chosen expansion point 
constrain the radius of convergence C. For the sharp cutoff and the expansion I, this has 



been analysed in [^] for = 1. It was found that t'sharp cannot be determined to an 
accuracy better than 8 ■ 10~^. Hence, the polynomial expansion I for a sharp cutoff does 
not converge beyond a certain level. 



However, the findings of ||32| do not imply that polynomial truncations are not trustwor- 
thy per se. To the contrary, the decisive difference between "good" or "bad" convergence 
properties stems essentially from an appropriate choice of the regularisation. When opti- 
mised, the regularisation implies a significant improvement for either expansion. In this 
light, the non-convergence of the sharp cutoff flow within expansion I is understood as a 
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deficiency of tlie sliarp cutoff regularisation. Tfiis picture is consistent witli tlie results of 
P^l , wlio sliowed tliat expansion II leads to an improved convergence over expansion I — 
even for the sharp cutoff. 



V. BOUNDS ON CRITICAL EXPONENTS 

In this section, we discuss how the critical exponents computed in the previous sections, 
depend on the regularisation. This discussion is mandatory because observables computed 
from a truncated flow, are known to depend spuriously on the regularisation. It is decisive to 
understand the range over which uerg may vary as a function of the IR regulator. The origin 
of the spurious RS dependence is easily understood. The regulator, while regulating the 
flow, also modifies the coupling amongst all vertex functions of the theory. These regulator 
induced contributions are of no relevance for the integrated full flow, but they do matter 
for approximated flows, like in the present case. It is argued that the smallest value for the 
exponent z/^j^^ is obtained for the optimised regulator -Ropt- Prior to this, we recall the results 
obtained previously in the literature, where, by a number of groups [p!3| , p2| , p0| , p2| , p4| , p!0 



critical exponents have been computed based on ( |2.1| ) for different regulators to leading 
order in the derivative expansion. For all N, all previously published results obey 

I^sharp > Z^ERG ^ '^min • (5.1) 

The regulators studied in the literature cover the sharp cutoff and a variety of smooth cutoff 
(exponential, power-law), and classes of regulators interpolating between the sharp cutoff 
and specific smooth cutoffs. Most results have been published for the Ising universality 
class = 1. For any > 0, the smallest value z/min obtained in the literature is larger 
than the value Uopt- i^min > t'opt- For a detailed comparison of critical exponents to leading 
and subleading order in the derivative expansion, and a comparison to results from other 



methods and experiment, we refer to Ref. |10 . 



A. Upper boundary 

Now, we turn to a general discussion on the scheme dependence of u. At first sight, 
( ^.1|) suggests that the possible range for u is bounded from above and from below. Let 
us assess the two boundaries. We begin by showing that the inequality t'sharp > '^erg does 
not hold for generic regulator. Indeed, the upper boundary in ( |5.1|) can be overcome by 
choosing regulators which lead to a worse convergence than the sharp cutoff. To see this 
more explicitly, consider a class of regulators discussed in Appendix |B[ It is given by 

Ra{q') = a{e-q')e{e-q'), (5.2) 

and is a variant of the optimised regulator ( p.ll| ), to which it reduces for a = 1. For a —>■ oo, 
it corresponds to the sharp cut off. The regulator leads to an effective radius of convergence 
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Ca = a ioT a < 1, and Ca = 2^ for a > 1. Since ^ 1 for a ^ 1, we expect that the 
corresponding critical exponents will become large. Our results for Ra in ( p.2| ) are given in 
Fig. 12 (for more details, see Appendix |B| and Tab. 3). We find that Va is a monotonously 
increasing function for decreasing a < 1. In particular, for small a we find indeed that 
> '^sharp- Heucc, this rcsult confirms the above picture: Flows with a poor radius of 
convergence lead to large numerical values for u. 
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Table 2: Critical exponents u (Ising universality class) for the flows and various 7. 



1.05 
1 

0.95 

0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 



~i r 



large- 



ERG 



Ra '■ 


• 


R^ '. 





Rc '■ 


opt 



sharp 



opt 



■ phys 



0.2 0.4 0.6 0.8 1 

X 



Figure 12: Ising universality class. The criti- 
cal exponent for various classes of regulators. 
For display purposes, we use x = |(7 — |) 
for R^, X = for Ra and x = c for Rc- 
Boundaries: The full line (opt) corresponds to 
Rc and denotes the lower boundary, the upper 
full line (large-A^) denotes the upper bound- 
ary. For comparison: The dashed line (sharp) 
indicates the sharp cut off value, and the thick 
full line (phys) the physical value. 



Next, we turn to a class of regulators aimed at probing the maximal spread of fiows. 
We recall that, to leading order in the derivative expansion, all regulator-dependence is 
contained in the functions i{uj). In Ref. it has been shown that the function i{uj) 
decays at most as if the regulator, obeying the basic constraints ( p^ - (|2.4|) , is not a 
strongly oscillating function of momenta. Hence, regulators leading to a function 



[UJ 



'1+UJ 



|1~7 



(5.3) 



define for 7 = 2 a boundary in the space of regulators. This is the case for -Ropt- No regulator 
can be found such that 7 > 2 [24]. The proportionality constant in (|5.3| ) is irrelevant because 
it can be scaled away into the fields in ( |2.7|) . The second boundary is set by a mass term 
regulator -Rmass = k"^'- For a mass term, the corresponding fiow is a Callan-Symanzik fiow 
which is not a Wilsonian fiow in the strict sense. This comes about because the condition 



17 



( p.3| ) does no longer hold true for -Rmass in the limit ^ oo. In higher dimensions, this may 
lead to an insufficiency in the integrating-out of large momentum modes. Inserting Rmass 
into in d = 3 dimensions leads to ( |5.3| ) with 1 = \- This sets the second boundary. 



The regulator is explicitly known for the cases 7 = 2, |, 1 and \, and corresponds, re- 
spectively, to the optimised regulator, the quartic regulator -Rquart = k'^ I'fi the sharp cutoff 
and a mass term regulator -Rmass = k'^- For all other values of 7 G [\, 2], R^ can be recovered 
explicitly from l^ioj) p4|. In the present case, we only need to know that such regulators 
exist. 

We have computed the critical index u-^ for l-^iuj) with 7 G [5, 2], and our results are 
given in Tab. 2. In Fig. 12, the results are denoted by open circles for a; = §(7 — ^). As a 
result, the function z/^ increases for decreasing 7, z/-y > z/opt. In particular, once 7 < 1, the 
results obey > z/sharp- When 7 — ^, the eigenvalues at criticality approach their large- 
values 1/ — > 1 and uJn —>■ 2n — 1 for any A^. Note that the large- A^ limit is exact in that 
it is independent of the regularisation |]TU[. The large numerical value for z/^=i/2 is due to 
the deficiencies of the Callan-Symanzik flow. We conjecture that the large- A^ limit i^i^^^^j^ 
corresponds, for any A^, to an upper boundary for any regulator 

J^^un < ^, M = 1 • (5.4) 

ERG — large N V / 



B. Lower boundary 

Next, we assess the lower boundary. It would be important to know whether the opti- 
mised regulator leads to the smallest attainable value for u in the present approximation. 
We are not aware of a general proof for this statement. However, strong evidence is provided 
by studying alterations of the optimised regulator. We have done so for various classes of 
regulators, three of which are discussed here more explicitly. For more details, we defer to 
the Appendices and 0. In Appendix |B|, we employ variants of the optimised regulator, 
given by the class Ra of (|5l^), and by the class Rb defined as 

Rbiq'') = (k^ - q') 0ik^ ~ g') - \k^) + {hk" + (1 - 2h)i) e{\k^ - q^) . (5.5) 



for < 6 < 00. This class contains the optimised regulator ( p.ll| ) for 6=1. The limit 



h ^ 00 corresponds to a variant of the standard sharp cutoff. Another new class of regulators 
Rc is studied in Appendix y, where we allow for an additional fc-dependence on a mass scale 
within the regulator. We use the class 

Here, c is a free parameter and m| = t/^' (0 = 0) is the mass term at vanishing field. For 
c = 0, it turns into the optimised regulator ( |2.11| ). Due to the implicit scale dependence of 
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TUk on k, the corresponding flow equations are substantially different from the usual one. 
Most notably, they contain terms proportional to the flow of m^. 



The classes R^, Ra, Rb and Rc probe "orthogonal" directions in the space of regulators. 
R^ is sensitive to the analyticity structure of the flow, Ra and R^ are sensitive to alterations 
of the function r{y) in the low momentum regime, and Rc-, while keeping the shape of the 
regulator r{y) flxed, alters the implicit modiflcations due to an additional running mass 
term. As such, i?^, Ri, and Rc can be seen as variants of the optimised regulator. The class 
i?-y covers the largest domain of qualitatively different flows . 




Figure 13: The spread of the critical expo- 
nent V for various N to leading order in the 
derivative expansion. The upper bound is set 
by the large- limit (Callan-Symanzik flow). 
The sharp cut-off results (Wegner-Houghton 
flow) are given for comparison. 



N 



We have computed the critical exponent u for the Ising universality class for all these 
classes of regulators. A part of our results for = 1 is displayed in Fig. 12. Similar re- 
sults are found for all A^. The classical (mean fleld) value for v is z/mf = \. The results 
for Ra in ( |5.2| ) are displayed in Fig. 12 by full circles, with x = Some numerical val- 
ues are collected in Tab. 3 (cf. Appendix P). It is found that all critical indices Pa obey 
^largoAT ^ ^ Kpt > ^mi- This proves that alterations of the regulator in the low momen- 
tum region do not lead to values for v smaller than Vopt- An analogous result is found for all 
regulators Rb in (|5.5| ), where z^i^^^g^^ > t'b > J-'opt- Some numerical values are given in Tab. 4 
(Appendix IB]). In Fig. 12, our results for Rb are represented by open circles, and those for 
Rc by a dashed line, with x = c. All regulators Rc from Eq. ( ^.6|) lead to the same critical 
exponent as -Ropt in Eq- ( p.llQ , I'c = ^opt (see Appendix y). 

In Fig. 13, we discuss the spread of u^^^^^ to leading order in the derivative expansion for 
all > 0. The spread i^,^^ / i^opt — 1 is a A^-dependent quantity. For A^ = 1, the spread is 
about 0.54, and hence quite large. For comparison, the relative width with respect to the 
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sharp cut off t'sharp/t'opt — 1 is roughly 0.06 and significantly smaller. With increasing A^, 
the spread vanishes as ~ This follows trivially from the fact that the ERG flow, to 

leading order in the derivative expansion, becomes exact in the large- limit [|T^. 

In summary, the critical exponent z/, as a function of the infrared regularisation, is 
bounded. The upper boundary is realised for flows with a mass term regulator, e.g. Callan- 
Symanzik flows. The lower boundary is given by 

z^ERG > ^opt (5.7) 

to leading order in the derivative expansion. Hence, z/opt appears indeed to be the smallest 
value attainable within the present approximation. The inequality ( ^.7|) provides a quan- 
titative basis for the optimisation procedure which has lead to Ropt in the first place. For 
the observable u, we have equally shown that the optimised regulator corresponds, at least, 
to a "local minimum" in the space of all regularisations. Furthermore, we have established 
flat directions in the space of regulators. Based on the conceptual reasons which have lead 
to Ropt expect that Uopt even corresponds to the global minimum. At present, 

however, we have no regulator-independent proof for this conjecture. 



VI. DISCUSSION AND CONCLUSIONS 

We studied in detail 0{N) symmetric scalar theories at criticality, using the ERG method 
to leading order in the derivative expansion. This included a complete investigation of the 
Gaussian flxed point in c? > 2 for arbitrary regulator, and the computation of universal 
critical exponents and subleading corrections at the Wilson-Fisher flxed point for d = 3. 
Furthermore, we studied the spurious scheme dependence for the critical exponent i> at the 
Wilson-Fisher flxed point in three dimensions. One of the main new results is that the 
leading critical exponent, as a function of the IR regulator, is bounded from above and from 
below as 

l^opt < Z^ERG ^ '^la.gc^ • (6-1) 

This result has been achieved by studying the maximal domain of ERG flows in the present 
approximation, ranging from Callan-Symanzik flows to optimised flows and variants thereof. 
The qualitative result - the existence of a non-trivial Wilson-Fisher flxed point - is very 
stable, although the spread of values for z/ is fairly large (see Fig. 13). Supposedly, this is a 
consequence of a small anomalous dimension t], which constrains higher order corrections. 
The spread would shrink to zero only to sufficiently high order in the derivative expansion 
or in the large- limit. 

The important quantitative question is: Which value for u could be considered as a good 
approximation to the physical theory? In view of the regulator dependence, a prediction 
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solely based on ( |6.1[ ) is of little use. Our answer to this problem is entirely based upon the 
structure of the ERG flow. We proposed to use specific regulators which lead to more stable 
ERG flows in the space of all action functionals. The numerical determination of critical 
exponents and subleading corrections to scaling as given in Tab. 1, is based on an optimised 
flow. We expect that the results should be in the vicinity of the physical theory. In the 
present approximation, the results for Uopt are indeed closest to the physical ones 

I^phys < ^opt ■ (6.2) 

The understanding of the spurious scheme dependence reduced the ambiguity in i/ to a small 
range about z/opt. Typically, the results from optimised flows other than i?opt are close to 
the values achieved by -Ropt, and hence close to the lower boundary of (|0| ). In this light, 
optimised flows are solutions to a "minimum sensitivity condition" in the space of all IR 
regularisations 0. 

Next we turn to the Callan-Symanzik flow, which is the flow with a mass term regulator 
-Rmass = k"^- We argued that it defines the upper boundary of values for z/. It is quite 
remarkable that this flow, for any A^, leads to the same eigenvalues at criticality given by 
the large-iV result. Here, this result has been achieved numerically. It would be helpful 
if it could also be understood analytically. The large numerical value for u reflects the 
poor convergence properties of a Callan-Symanzik flow, essentially due to deficiencies in the 
integrating-out of large momentum modes. 

Now we discuss our results concerning polynomial approximations. The reliability of 
this additional truncation is guaranteed if the approximation convergences reasonably fast. 
Here, we have established that optimised flows converge very rapidly within the local po- 
tential approximation. The efficiency is remarkable: a simple approximation with only six 
independent operators — say, the running v.e.v. and five running vertex functions up to 
(0"0a)^ — reproduces the physical result for the exponent z/ at the percent level. A bet- 
ter agreement cannot be expected, given that anomalous dimensions of the order of a few 
percent have been neglected. These findings are in contrast to earlier computations based 
on the sharp cut-off, where the polynomial approximation has lead to spurious fixed point 
solutions, even to high order in the approximation Jl^. Hence, the efficiency of the for- 
malism not only depends on the choice for the degrees of freedom and the truncation, but 
additionally, and strongly, on the IR regulator. We conclude that polynomial approxima- 
tions are reliable for all technical purposes, and even to low orders, if they are backed-up 
by appropriate regulators. These considerations should be useful in more complex theories 
whose algebraic complexity requires polynomial approximations, e.g. quantum gravity. 

It would be interesting to apply the present ideas to theories like QCD, where the 
propagating modes are strongly modified in the low momentum regime due to confinement 
36| , |37[] . Then, an optimised regulator is found by requiring that the regularised inverse prop- 



agator is again flat, i. e. momentum-independent for small momenta. Interesting choices 
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are R = [hP' — X) 6{k'^ — X) and X = [0 = 0o] and variants thereof. Here 0o denotes a 
non-propagating background field. This conjecture is supported by first results. 
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A. GAUSSIAN FIXED POINT 

In this appendix, we discuss the Gaussian fixed point of the flow equation ( |2.13D in > 2 
dimensions and for arbitrary regulator. The Gaussian fixed point corresponds to the specific 
solution u^,{p) = const. All higher derivatives of the potential vanish, u^\p) = 0. From the 
flow equation, we deduce that 



Jo l + riy) 



d Jo 1 + r{y) 



For the optimised regulator, we find = ANvd/d"^. More information can be extracted 
by studying small perturbations Su^"^^ around the m-th derivative of the scaling solution 
Ui,{pY'^\ The eigenperturbations obey the differential equation 

dt Su^""^ = u Su^"''^ (A.2) 

with eigenvalues u. Expanding the flow equation to leading order in 6u, the eigenvalues 
obey 

0=[uj + d~{d- 2)m] Su^"^^ + 2A(y + m)-{d~ 2)p Su^""^^^ + 2pA Su^"'^^^ . (A.3) 
Here, the scheme- dependent coefficient A is given by 



A = 2vJ dyy2~^—-^f^ (A.4) 



[l + r{y)Y 

and < A < oo. For the optimised regulator, Aopt = 4:Vd/d. Introducing new variables 
X = {d — 2)p/{2A) and /(x) = Su{p), the differential equation ( |A.3| ) transforms into the 
(generahsed) Laguerre differential equation 



U-2 



/(™) {x)+(^^ + m- /(™+^) (x) + (x) . (A.5) 



We consider only polynomial solutions to Eq. ( |A.5| ). The requirement that solutions to 
Eq. ( [A.5| ) are bounded by polynomials fixes the possible eigenvalues as 



22 



00 



[d - 2){n + m) - d 



(A.6) 



for non- negative integers n and m. Apart from an irrelevant normalisation constant, the 
77,-tli eigensolution to Eq. (|A.5|) are given by the (generalised) Laguerre polynomials 



2Ap 
d-2 



(A.7) 



Eq. ( |A.7| ) is the most general eigensolution at the Gaussian fixed point in d > 2 dimensions 
with eigenvalues given by Eq. ( |A.6| ). The result holds for arbitrary regulator function. The 
scheme dependence enters only the argument of the Laguerre polynomials in Eq. ( |A.7| ). 
It is interesting to note that the eigenvalues are independent of the regularisation scheme. 
Furthermore, the rescaled differential equation (|A.5| ) is also independent of the regulator. 



In d = 3 dimensions, the relevant and marginal operators are Lq 



N/2-1 



and L 



N/2-1 



jN/2-1 
^1 1 



r N/2-1 
^2 



with eigenvalues uj = —3, —2, —1 and for m 



0, or Lq^'^, Li^"^ and L^^^ with 



eigenvalues ui = —2, —1 and for m = 1. In d = 4 dimensions, the relevant and marginal 



operators are Lq 



N/2-1 



Li^'^ ^ and L^^^ ^ with eigenvalues oo = —4, —2 and for m = 0, or 



Lq^'^ and L^^'^ with eigenvalues u = —2 and for m = 1. For m > all eigenoperators 
are marginal or irrelevant. 



For 2m + = 1 (2m + A^ = 3), the solution Eq. ( A.7 ) can be rewritten in terms of 
Hermite polynomials of even (odd) degree. 



m 



m 



l-N 



3-N 




(A.8) 
(A.9) 



Let us finally mention that some of these solutions have been given earlier in the literature 
for the case of a sharp cutoff regulator (with Agharp = 1)^ for A^ = 1 and m = 1, Eq. ( |A.9| ) 
has been given in Ref. |T3], and Eq. (|A.7|) for m = 1 has been given in Ref. 



B. VARIANTS OF THE OPTIMISED CUTOFF 

In this appendix, we discuss variants of the optimised regulator (|2.11|) . The aim is to 
probe whether certain alterations of the regulator may lead to lower values for the critical 
exponent v. Here, the properties of the regularisation are changed in the low momentum 
region by modifying the function riy). In the following appendix, we discuss modifications 
of r{y) through the introduction of additional (theory-dependent) k dependent parameters. 
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Figure 14: The function y{l+ra) for the class 
of regulators ( [B.l| ) . All lines for different a co- 
incide for large momenta y G [1, oo] (full line), 
but differ for small momenta y £ [0, 1[ (dashed 
lines). The dashed lines clock-wise from the 
bottom: a = ^, 1, 2 and oo. The short dashed 
line corresponds to a = (no regulator). 



1. Definition 



In this appendix, we discuss two variants of the optimised regulator ( p.ll| ). First, we 
consider the class of regulators given by 



(B.l) 



These regulators have a compact support. They vanish for all > k"^. In the infrared, 
they have a mass-like limit Ra{q^ — > 0) = ak"^ for all < a < oo. The limit a — > cxo 
corresponds to the sharp cut-off case. For a = 0, the regulator is removed completely. 
For a = 1, the regulator ( [B.l| ) reduces to the optimised regulator (|2.11| ). As a function of 
a, the regulators differ only in the momentum regime y G [0,1[ , where y = q^/k"^. The 



dimensionless functions corresponding to ( p.l|) are given by 



r„(y) = a (i- 1)^^(1 



y) 



(B.2) 



In Fig. 1, we have displayed the function y{l + r^) for various cases. The full line corre- 
sponds to the range y G [l,c>o], for all regulators ( |B.1| ). The dashed lines, clock-wise from 
the bottom, correspond to a = ^,1,2 and oo. The gaps associated to ( B.l ) are given by 
Ca = 2^±i for a > 1, and Ca = a for ^ < a < 1. They are obtained from the normalised 
analogue of ( p.l|) , chosen such that 



1. 



Second, we consider another variant of the optimised regulator, where the properties of 
the regularisation are changed only in the low momentum region by modifying the function 
r{y). Consider the class of regulators given by 



i?fe(g2) = (^2 _ ^2) Q^^2 _ ^2) q^^2 _ 1^2) ^ ^^^2 ^ _ ^^^^2^ ^^1^2 _ ^2 



(B.3) 
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These regulators have a compact support. They vanish for all > k"^. In the infrared, they 
have a mass-like limit Rb{q^ —>■ 0) = bk"^ for all < 6 < oo. At first sight, it may seem that 
( [B.3D is not a viable regulator for 6 = 0, because Rb=o vanishes in the infrared limit (no gap). 
However, for 6 = the function dtRb=o vanishes identically for all q"^ < Hence, ( p.3|) 
provides a gap because Rb=o{^k'^) = \k'^ > 0. For 6=1, the regulator ( [B.3|) reduces to the 
optimised regulator (|2.11|) . As a function of b, the regulators differ only in the momentum 
regime y G [0, |[ , where y = q^ /k"^. The dimensionless functions Vb corresponding to ( p.3| ) 
are given by 



n{y) 



■-y ' 



26) Q{\ - y) . 



(B.4) 



In Fig. 15, we have displayed the function y{l+rb) for various cases. The full line corresponds 
to the range y G [^,C)o], for all regulators ( p.3|) . The dashed lines, clock- wise from the 
bottom, correspond to 6 = 0, i, 1, 2 and oo. By construction, the regulator is normalised as 
rft(^) = 1. The associated gaps are given by = 1 for 6 > 1 and 6 = (see below), and 
Cb = b for < 6 < 1. For comparison, we have also given the curve for the standard sharp 
cutoff (dotted line). The corresponding gap is Cgharp = ^• 




Figure 15: The function y{l+rb) for the class 
of regulators ( |B.3| ). All lines for different b co- 
incide for large momenta ?/ E [^i oo] (full line), 
but differ for small momenta y G [0, ^] (dashed 
lines). The dashed lines clock-wise from the 
bottom: b = 0,^,1,2 and co. The standard 
sharp cutoff regulator (dotted line) is given for 
comparison. 



The regulators Ra and Rb have a similar low momentum limit for 0, e.g. Ra{q^ 
0) = Rb=a{q^ — >■ 0) = ak'^. The crucial difference between them concerns the intermediate 
momentum regime q^ ~ fc^. Here, the regulator Rb leads by construction to a plateau for 
y{l + Vb), which is absent for Ra- 

2. Flows 



In order to employ (^ 
fiows £(u!) 



for the computation of critical exponents in d = 3, the associated 
u) have to be computed. For the function ia{^) and for a < 1, we find in 
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d = 3 (similar expressions are found for any d) 



2a 



la + iu 



1 



arctan 



a 



/I - a 
a + Lu 



(B.5) 



The region a > 1 is obtained by analytical continuation: 

2a 



1 



la + 00 



a 



arctanhi 



a + u 



(B.6) 



These flows have the following structure. They have poles on the negative u-axis at u = —1. 
For a = oo, the function decays only logarithmically for asymptotically large uj. In the 
limiting cases a = 1 and oo, we find 



1(1+^) 



ln(l + uj) + const. 



(B.7) 
(B.8) 



For a = 1 , in ( [B.7| ) , both expressions ( |B.5| ) and ( [B.6D have the same limit discussed earlier in 
P]. In the limit a — > oo, the resulting regulator is equivalent to the standard sharp cut-off. 



In order to employ the regulator Rb from ( |B.3| ) for the computation of critical exponents 
in d = 3, the associated flows i{uj) h{^^) have to be computed. In full analogy to the 
preceding computation, we find in d = 3 (similar expressions are found for any d) for the 
function ib{uj) and for 6 < 1 



[l-2-^/^){l + 



UJ 



1 



Ib + u 



V2{1 - b) 

The region 6 > 1 is obtained by analytical continuation: 

h 



arctan 



1 



(B.9) 



4(cu) = f(l-2-3/2)(l + ^)-i + 



72(1 - h) 



-arctanhi 



'6-1 



(B.IO) 



These flows have the following structure. They have poles on the negative tu-axis at = — 1 
due to the first term on the r.h.s. in ( |B.9| ) and ( [B.lOp . For large uj and h < oo, both 
expressions decay as uj~^. For b = oo, the decay is only logarithmic. Let us consider the 
three limiting cases a = 0, 1 and oo. We find 



ib=oH = 1(1 - 2~'/'){i + ujy' 

h=iiuj) 



3* 

1(1 



UJ 



-'6=00 



.UJ 



[1 - 2-3/2)(i + _ _i_ in(l + uj) + const 



(B.ll) 
(B.12) 
(B.13) 



For 6 = 0, the momentum regime < ifc^ does not contribute to the flow and the effective 



gap for 6 = is Co = 1. This is seen directly from 



and (B.ll): the first term of 
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^ stems from the momentum interval ^k^ < < k"^, which is the only term surviving 
For 6 = 1, in { B.V2 ), both terms of ( [B.9| ) combine to the known result discussed 
51. Finally, we turn to the limit b oo. The resulting regulator is similar to the 



in (|KTT]) . 

earlier in 

standard sharp cut-off, with, however, an important difference. For the sharp cutoff, the 
function y{l + ri,) has no plateau in the momentum regime ^k"^ < < k"^ (see Fig. 15), 
which leads to 4harp(t^) = — ln(l -l-a;). In (|B.8|) , the sharp-cutoff-like logarithmic term is 
clearly seen, and is due to the momentum integration with y G [0, However, a decisive 
difference is the additional term in (p.8| ). Notice also that the constant in ( p.l3|) is actually 
infinite, but field independent. Hence, it is irrelevant for a computation of critical exponents 
(only the functions d^jtbi'^) are needed). 



3. Results 



For we have computed the critical exponent u for the Ising universality class using 
the flow equation ( |2.7|) with (|B.5| ), (|B.6| ) (see Tab. 3). We confirm that i'a=oo = t'sharp and 
that Va=i = Vopf For all a > 1 (a < 1), z/^ is a monotonically increasing (decreasing) 
function with increasing a, hence z/q > z/opt- Notice that the smallest value for v is obtained 
for the largest value for the gap parameter. 



a 


10-2 


10-1 


1 

2 


1 


2 


10 


102 


10^ 


10^ 


oo 


Ca 


10-2 


10-1 


1 

2 


1 


2 
3 


10 
19 


100 
199 


1000 
1999 


10000 
19999 


1 

2 




.776 


.677 


.652 


.650 


.651 


.665 


.683 


.688 


.689 


.690 



Table 3: Critical exponents u (Ising universality class) for the regulator Ra and various a. 



b 





1 

2 


1 


10 


100 


oo 




.6495 


.6518 


.6495 


.6594 


.6675 


.6699 



Table 4: Critical exponents v (Ising universality class) for the flows with and various h. 

For i?f,, we have also computed the critical exponent v for the Ising universality class 
using the flow equation (|2.71) with ([B.9[ ), (IB.lOl) . We find that z/b=o = i^6=i = t'opt- For 



6 < 1, we have z/f, > z/opt- However, for too small values of 6, 1 ^ 6 > 0, the regulator 
leads to a very small gap and the polynomial approximation does no longer converge to a 
definite result, which is an artifact of the regulator. For 6 > 1, z/;, is a monotonic function 
of b with z/fe=oo > J^b > t^opt- It is interesting to discuss the case 6 = oo in more detail. Here, 
h'b=oo = -6699 ■ ■ ■ which should be compared to z/gharp = -6895 ■ ■ ■ and to z^opt = -6495 ■ ■ ■. 
From the structure of the flow, it is clear that the difference i'sharp — ^b=oo has to be attributed 
to momenta with y G [^,1]- Hence, the "flattening" of the standard sharp cutoff reduces 
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the critical exponent by a few percent. On the technical level, this reduction is attributed 
to the non-logarithmic term in ( [B.13| ). However, the smallest value for v is obtained only 
in case the logarithmic term is absent, as it happens in both (|B.11|) and (|B.12|). 



C. CUTOFFS WITH INTRINSIC SCALING 

In this appendix, we study regulators with additional intrinsic scale-dependent parame- 
ters. 



1. Definition 

Up to now, we have considered IR regulators R{q^) which depend on momenta only 
through the combination /k"^, cf. Eq. (p75|). In Eq. ( |2.5| ), the essential IR cutoff is provided 
by the function r{q^ /k'^), which cuts off the momentum scale q^ in a way which is independent 
of the particular theory studied. The k dependence of R{q^) is the trivial k dependence hnked 
to the dimensionality of R. The situation is different once further fc-dependent functions 
are introduced into the regulator. Typically, this is done by replacing 

R{q')^R{q\{Zu,ml...}). (C.l) 

Here, the set {Z^, m^, . . .} denotes scale-dependent parameters of the specific theory studied, 
like the wave-function renormalisation Zk or mass parameters m^. It is expected that the 
substitution ( |C1[ ) leads to an improved convergence and stability of the flow. A well-known 
example for ( |(J.1D is given by R{q^) — ^ Z^Ri^q^) , which is often used beyond the leading 
order in a derivative expansion. Here, the introduction of Zi^ in the regulator simplifies 
the study of scaling solutions. In the context of non-Abelian gauge theories, regulators like 
Eq. (|C.1|) have been used in |^ based on the replacement q^ ^^^''[^o] in the regulator, 
and hence 

R{q') i?[f f ] . (C.2) 

Again, this type of regulator has been motivated to stabilise the flow and to encompass 
possible poles in the flow due to mass terms for the gluonic fields . Notice that T). in 



Eq. (|C.2D cannot be the full field-dependent functional, because elsewise the flow equation 
would no longer be the correct one. Instead, one has to evaluate it for some fixed background 
field = 00 (hence the bar on f [f"*). In the present case, and to leading order in the derivative 
expansion, we have 



ff [0o] = g^ + f/n0o). (C.3) 

Here, f/"(0o) corresponds to a scale-dependent effective mass term. Within the non-convex 
part of the potential, or in a phase with spontaneous symmetry breaking, we have m| = 
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U"{(f)o) < 0. Hence the regulator Eq. ( p.2| ) is a special case of Eq. ( |C.1D . In the remainder 



the optimised regulator ( |2.12| ) is used to define a class of regulators of the form ( |(J.1|) , namely 

Hence, Rc{<f') = Ropti.'f' + cm\). Effectively, this corresponds to the replacement k"^ — >• 
^eff(^) = + cm\. In terms of r{y) defined in Eq. ( p.5|) , the regulator is given as 

rc{y) = {^-l)e{l-y-cu), (C.5) 

and the dimensionless mass parameter is 

u = ml/k\ (C.6) 

Below, we use ml = U"{(f)o = 0). The regulator ( |(J.4| ) can be seen as a 'sliding' cutoff, 
because at any scale k, only the fc-dependent momentum interval y G [0, l—cCOk] contributes 
to the flow. For c = 0, the regulator reduces to while for c = 1 it turns into a regulator 



of the form (|C.2|) with ( |C.3|) . Due to the additional dependence on m^/A;^, the structure of 



the flow equation is now different. 



2. Flows 

Let us compute the flow for the regulator Eq. ( |C.4| ). Inserting Eq. ( |C.4| ) into Eq. 
and after some straightforward algebra, we flnd 

f.M-id-c^j-'/^ifl^ (C.7) 

i — C UJ -\- UJ 

Notice that the function Eq. ( |(J.7| ) depends now on u and on the flow of u. This is generic 
to regulators of the form (|C.1D or (|C.2|) , because the implicit scale dependence of ml leads 
to an additional term ~ dtml in the flow equation. The flow equation ( |2.7|) with Eq. ( |C.7| ) 
becomes 



dtu + du-{d- 2)pu' = {N- 1)(1 - cCu) 



d/2 l-CU-§dtUJ 

1 — cu + u' 



^ ^ 1- CLU + u' + 2pu" ^ ' 

Here, in order to simplify the expressions, we have rescaled the irrelevant numerical factor 
into the flelds and the potential. In order to flnd an explicit form of the flow, the running 
mass term needs to be specifled. We chose u = u'{p = 0). This choice is motivated by the 
fact that the function u', on the flxed point, reaches its most negative value at vanishing 
fleld. For c = 0, the original flow may run into a pole at u' = —1. The present choice shifts 
the pole to u'{p) = — H-cm'(O). Since u'{p) —u'{0) > for a scaling solution, the right-hand 
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side of ( p.8| ) has no longer an explicit pole for c = 1. This has been the motivation for the 
structure of the regulator used in However, as we shall see, the full flow still has an 
implicit pole due to the flow of uj in ( |(J.8D . In terms of 



Ao = u'ip = 0) , Ai = u"ip = 0) 



(C.9) 



and after inserting Eq. (|C.7|) into Eq. (|C.8|) and collecting terms proportional to the flow of 
Ao, we find 



dtXc 



2Ao + Ai(iV + 2)(1 - c Ao)^+'^/^[l + Ao(l - c)] 
' 1 - icAi(iV + 2)(1 - c Ao)'^/2[l + Ao(l - c)]- 



-2 



(C.IO) 



For the quartic coupling, and suppressing terms proportional to u"'{0), we find 



dtXi 



-Ai Ai 



2{N + 8) (1 -cAq 
[l + (l-c)Ao]' 



X 



^ c Ao + ^(iV + 2)(l-cAo)^+'^/^[l + (l-c) Ao]- 
° 2 1 _ c^(iV + 2) (1 - cAo)'/'[l + (1 - c) Ao]-' 



(c.ii; 



The structure of the flows (|C.10|) and (|C.11|) is easily understood. The denominator in 
( p.lO|) stems from a resummation of the back coupling of dtXo- The denominator becomes 
trivial for c = 0. The numerator contains the usual scaling term and a modified threshold 
behaviour, which depends now on c. A similar structure appears for the flow of the quartic 
coupling. Simpler forms of the flow are obtained for c = (no back-coupling of a running 
mass term) or c = 1. For c = 0, we have 



dtXo = -2Ao + Ai 
dtXi-- 



N + 2 



yl + XoY 
(4 - d)Xi + Xj 



,2{N + 8) 



(1 + Xof 
while for c = 1 , the result is 

2Ao + Ai(iV + 2)(l-Ao)i+'^/2 



(C.12) 
(C.13) 



dtXo 



dtXi = -Ai(l + Ao 



iAi(iV + 2)(l-Ao)'^/2 ' 

4 - - Ai 2(iV + 8) (1 - Ao)'^/^ 



l-iAi(iV + 2)(l-Ao)^/2 



(C.14) 
(C.15) 



Hence, ( |C.12|) and ( |C.13| ) have a putative pole at Aq = — 1. In turn, ( p. 14]) and (|C.15|) have 
a putative pole at {N + 2)Ai = 2(1 — Ao)"'^''^. The putative pole is absent for N = —2. 
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3. Results 



From the explicit form of the function (|C.7|) , and without having made yet the exphcit 



choice (|C.9|) for the mass term ml, we conclude that the non- universal fixed point solution 



dtu = of (pl7| ) for either i?opt from (|2.11| ), or the class of regulators Rc from (|C.4|) are 



related by a simple rescaling of the fields. The reason is the following: on the fixed point 
dfU = 0, we also have dt{m\/k'^) = 0. Hence, the functions Eq. ( |C.7| ) do no longer depend 
on the flow of the mass term. Hence, Eq. ( p.7| ) becomes i{uj) = |C"^/^+^/(C + cu), with 
C = 1 — cuj = const on a flxed point. This function is equivalent to a flow derived from 
Ropt- Hence, it suffices to rescale u — > u/C^^"^ and p — > p/C^^^^^ in order to transform the 
flxed point solution for arbitrary c onto the flxed point equation for c = 0. 

Next, we check the c dependence of critical exponents. A priori, the critical exponents 
are sensitive to the flow in the vicinity of the flxed point, and hence to the additional terms 
dtXf) in the flow equation. We have solved the flow (|C.8|) with ( |C.1CI|) numerically, within 



the polynomial approximation deflned in Eq. (|3.1| ) up to n^rmc = 20. The eigenvalues at 
criticality are found to be independent on the parameter c. The results correspond to the 
lower full line in Fig. 12 (and x = c). Furthermore, we found that the critical exponents for 
different c agree for every single order in the truncation. Hence, for = 1, the exponent 
i^trunc as obtained from ( |U.8| ) with ( |(J.1(J| ) are given by the line I in Fig. 10. 

In summary, the introduction of scale-dependent parameters into the regulator has lead 
to a signiflcant change of the flow equations and their pole structure. Hence, the approach 
to a flxed point solution, and stability properties of the flow are quite different. Here, we 
studied the replacement Ropt{,<f) — ^ Ropti'f + cm\). Universal quantities are independent 
of c, because the modifled regulator is linked to the optimised one through the replacement 
k"^ — > k'^Q{k) = k"^ + cm^ in the flow equation. This should be irrelevant for universal 
observables at a flxed point, as has been conflrmed explicitly. For the same reason, the 
entire class of regulators Ra{q^ + cm\) = aRoptiq"^ + cm\) leads to c-independent critical 
exponents. Still, the flows are completely different for different a. Therefore, c can be 
used as a free parameter to stabilise the flow, without affecting the physical result. If the 
substitution R{q^) R{q^ +cm\) cannot be rephrased as a redeflnition of the infrared scale, 
it is expected that also universal observables no longer remain insensitive to free parameters 
like c. This case should be studied separately. 
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